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Abstract 

,S^ . We propose a formalisni for extracting molecular Hamiltonian structure from inversion 

rp-i' of time-dependent fluorescence intensity data. The proposed method requires a minimum of 

a priori knowledge about the system and allows for extracting a complete set of information 
about the Hamiltonian for a pair of molecular electronic surfaces. 
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1 Introduction 



O 

O 

c/2 I A long standing objective is the extraction of molecular Hamiltonian information from laboratory 

^^' data. The traditional approaches to this problem attempt to make use of time-independent 

O ■ (spectroscopic and scattering) data [|l|, |2|, ^. Another approach aims to use ultrafast temporal 

data, with information on molecular potentials and dipole moments obtained for spatial regions 
sampled by evolving wave packets. Research in this direction has been especially intense during 
K* ■ the last few years [y, ^, g, 0, |g, |§] • This activity is inspired by recent progress in the technology 

of ultrafast laser pulses |jl^, ^], which makes possible observations of molecular dynamics with 
f^ [ increasingly higher spatial and temporal resolution. 

^1 ' Due to the difficulty of the Hamiltonian inversion problem, it is common to assume that 

some a priori knowledge of the system is available. For example, one technique [y] proposes 
f^ I to extract time-evolving wave functions and excited-state potentials using time-resolved and 

c/3 ■ frequency-resolved fluorescence data and knowledge of the ground-state potential, the transition 

frequencies, and the transition dipole moment. The inverse tracking method Q, proposed 
£^| for recovering the potential energy and dipole moment of a molecular electronic surface by 

r~| ' monitoring the temporal evolution of wave packets, explicitly assumes knowledge of the initial 

Qh' excited wave functions. Clearly, such assumptions impair self-consistency and at least partially 

undermine the inversion objectives. Although the desire to simplify the inversion algorithm 
K^ \ by making a priori assumptions about what is known and unknown is understandable, it has 

^ • remained an open question about whether these assumptions are actually necessary. 

^ This letter addresses the latter point by proposing an inversion formalism that makes use of 

minimal a priori knowledge about the system. The formalism is designed to operate between 
two electronic surfaces, with electronic and vibrational transitions driven by two fast laser pulses, 
which allows for extracting the potential energies and dipole moments for both surfaces as well 
as the electronic transition dipole moment. The extraction is based on the inversion of the 
time-dependent fluorescence intensity data obtained from the detection of spontaneous emission 
in transitions between the electronic surfaces. The proposed formalism lays the ground work 
for extracting a complete set of information about a pair of electronic surfaces in a closed way, 
with a minimum of a priori assumptions about the molecular Hamiltonian. This letter presents 



the conceptual foundation of this novel approach, and a detailed numerical algorithm with 
simulations will be presented elsewhere. 

2 The physical picture 

Consider the ground and excited electronic molecular surfaces with potential energies Vg{x) and 
Ve{x) and dipole moments Hg{x) and ^e{x)^ respectively. The dipole moment for the electronic 
transition between the two surfaces is M{x). For the sake of conceptual clarity, we consider a 
one-dimensional problem; the generalization for the multidimensional case is straightforward. 

The setup includes two time-dependent locked laser fields: eo(t) drives transitions between 
the two electronic surfaces (the carrier frequency of this laser will be typically in the visible or 
ultraviolet part of the spectrum), and ei(t) drives transitions between vibrational levels within 
each of the two surfaces (the carrier frequency of this laser will be typically in the infrared). 
The role of the driving fields is to excite the molecular wave packet and guide its motion on 
the surfaces. It is physically reasonable that the potentials and dipole moments may be reliably 
extracted only in the region sampled by the evolving wave packet. We assume that interactions 
with other electronic surfaces and incoherent processes (e.g., thermal excitation and collisional 
relaxation) are negligible. Ultrafast laser technology has made great advances recently, but 
preparation of the infrared pulse of a desired shape is still a challenging technical problem. We 
will consider the general situation, with two locked driving fields and five unknown functions 
(two potentials and three dipole moments), but taking ei = the problem is easily reduced to 
a simpler one, with only one driving field eo and three unknown functions (VJy, Vg and M). 

The Hamiltonian of the system in the Born-Oppenheimer, electric-dipole and rotating-wave 
approximations takes the form: 

H = Hg{x,p,t)agg + He{x,p,t)aee - M{x)€o{t){aeg + (Jge), (1) 

where aij = \i){j\ (with i,j = g,e) are transition-projection operators for the electronic states 
\g) and |e). Here, x and p are the canonical position and momentum for the vibrational degree 
of freedom, Hg and Hg are the vibrational Hamiltonians in the ground and excited electronic 
states, 

Hi{x,p,t) =T + Vi{x) - Hi{x)ei{t), i = e,g, (2) 

and T = p^ /2m is the kinetic energy of the vibrational motion. 

We assume that the initial state of the system is |^(0)) = \uQ)\g), where uq{x) is the 
vibrational ground state localized in the known harmonic part of the potential Vg{x). The 
state of the system at any time t will be of the form 

\^(t)) = \u{t))\g) + \v{t))\e), (3) 

with the normalization condition J dx (|n(x, t)p -|- |t;(x,t)p) = 1. The Schrodinger equation, 
ihdt\^{t)) = H\^{t)), then takes the form 

ihdtu{x,t) = -- — dlu{x,t) + \Vg{x) - fig{x)ei{t)]u{x,t) - M{x)eo{t)v{x,t), (4) 

ihdtv{x, t) = --—dlv{x, t) + \Ve{x) - ^ie{x)ei{t)]v{x, t) - M{x)eo{t)u{x, t), (5) 

2m 

with the initial conditions u{x,0) = uo{x), v{x,0) = 0. 



The radiation emitted spontaneously by the molecule via transitions between the excited and 
ground electronic surfaces contains information about the wave packet. This fact has been used 
to reconstruct unknown vibrational wave packets in the method of emission tomography [p^]. 
Our aim is different: we assume that the initial state of the system is known and want to extract 
the unknown potentials {Vg and Ve) and dipole moments (ng, fie, and M) from information 
contained in the time-dependent fluorescence. We choose the time-dependent intensity of the 
emitted radiation, I{t), as the observable. This intensity is I{t) = E^^\t)E^' {t) , where E^~^\t) 
is the negative-frequency part of the electric field operator of the emitted radiation. E^^' is 
proportional to Ma eg, so the measured quantity is 

(I(t))=Ac(^(t)|MVee|^(t)). (6) 

where k is a proportionality constant. 

3 Extraction of the Hamiltonian structure 

The physical picture above leads to the the following mathematical problem: extract the po- 
tentials and dipole moments from the measured intensity {I{t)), assuming that the initial state 
and the two driving fields are known (note that a number of advanced experimental techniques 



have been recently developed for characterization of ultrashort optical pulses [14, |l^, |l 



We start from the Heisenberg equation of motion, ihdi/dt = [I, H], to obtain 

{^{t)\[M\T]aee-eo{t)M^{aeg-age)\^{t)). (7) 



K dt 
Using form (H) of the wave function, we rewrite (1^) as an integral equation for M{x): 

dx[M\x)F{x, t) + M^{x)G{x, t)] = A^, (8) 

2k dt 

where 

F(x, t) = Im[t;(x, t)dlv*(x, t)], G(x, t) = eo(i)Im[n*(x, t)v{x, t)]. (9) 

2m 

In order to obtain equations for the other unknown functions (two potentials and two dipole 
moments), we consider the second time derivative of {I{t)). Then, using (y), we derive the 
following integral equation: 

h^ ^2 / j\ r 

— ^^ + ^M{t) = dx[Ve{x)-Vg{x)-ei{t)iJe{x)+ei{t)fig{x)]SM{x,t) 



+ / dx[Veix) - ei{t)fle{x)]RM{x,t), (10) 



where 



??2 

RM{x,t) = —Ke{v*{x,t)[dl,M'^{x)]v{x,t)] , (11) 

SM{x,t) = -2eo{t)M'^{x)Re[u*{x,t)v*{x,t)] , (12) 



fr r 

TmU) = —^dxv*{x,t)\M'^(x)dt-2dlM'^ix)dl + dtM'^ix)]v{x,t) 
A.ni'^ J 

+— eo(t)Re| I dxv*{x,t)[M^{x)dl + M'^{x)dlM{x) -2dlM^{x)]u{x,t) 
+2el{t) I dxM^{x) (|u(x,t)|2 - \u{x,t)\^) . (13) 



It is convenient to formally enumerate the unknown functions, 

/i(x) = Vg(x), f2{x) = Ve{x), h{x) = doHg{x), /^{x) = dofieix), (14) 

where do = 1 V/m, so all the functions fr{x) have the dimension of energy. Then the integral 
equation ( [Tol ) takes the form 

r ^ 

j dxY,Kr{x,t)U{x)=g{t), (15) 



where 

Ki{x,t) = -SM{x,t), K2{x,t) = RMix,t) + SMix,t), (16) 

K3{x,t) = -ei{t)Ki{x,t), Ki{x,t) = -ei{t)K2{x,t), (17) 

git) = ^^+TM{t), (18) 

and ei = dQ ei is the scaled (dimensionless) field. 

It is important to emphasize that in fact equations (^) and ( |T^ ) represent an infinite number 
(or, in practice, a large number) of equations corresponding to different times. We will use this 
fact in the regularization procedure below. Of course, equation ( [l5| ) is nonlinear because the 
wave function depends on the potentials and dipole moments. Similarly, a solution M of (^) 
depends on the wave function and thereby depends on other unknown functions. Consequently, 
the problem at hand, including the integral equations and the Schrodinger equation, is highly 
nonlinear. More importantly, the solution for such a system of integral equations is generally not 
unique and the problem is ill-posed (i.e., the solution is not stable against small changes of the 
data). These characteristics are common to virtually all inverse problems and arise because the 
data used for the inversion are inevitably incomplete. Consequently, we need to regularize the 
problem by imposing physically motivated constraints on the unknown functions. For example, 
we may use the fact that physically acceptable potentials and dipoles should be smooth functions 
and tend to zero asymptotically as a; ^ oo (in the case of the dipole, the atoms are assumed to 
separate as neutrals). By taking into account this information, some constraints are imposed on 
the solutions, singling out the functions with desirable physical properties. This regularization 
procedure will stabilize the solution. 

The regularized solution of equation (|l5|) is achieved by minimizing the functional 

dxf^{x). (19) 

r=l 



, r 4 "I 2 4 

J= f dt' fdxY, Kr{x, t')fr{x) - g{t') +^ar f 

•^0 l-J j.=l J r=l -^ 



Here, ar are standard regularization parameters which denote the tradeoff between reproducing 
the laboratory data and obtaining the solution with smooth and regular functions. The time 



integration in (|l9| ) has a simple physical meaning: the measured intensity brings in information 
about the potentials and dipoles at each instance of time and we want to use all the laboratory 
information which has been accumulated during the period from time zero until t. The choice of 
the functional (|l^ ) is not unique, and other forms of regularization may be considered as well. 
Taking the variation of the functional J^ with respect to the unknown functions frix) involves 
a subtlety related to the nonlinearity of the problem: the kernels K^ and the free term g depend 
on the wave function and on M{x) and thereby depend on fr{x). The practical (numerical) 
solution of any nonlinear problem includes some kind of linearization. Here, the point at which 
we make the linearization is determined in the regularization procedure. We choose to take the 
variation of the functional J in equation ([l9|) only with respect to the explicit dependence on 
frix). Then we obtain the set of regularized equations: 



where 



/ dx' ^ ^r{x, x', t)fr{x') + apfp{x) = 5p(x, t), 

J r=l 



M^r{x,x',t)= / dt'Kp{x,t')Kr{x',t'), 



9p{x,t)= / dt'Kp{x,t')g{t'). 
Jo 



(20) 

(21) 
(22) 



With p,r = 1, 2, 3, 4, we have the system of four integral equations with four unknown functions 
(two potentials and two moments). 

Now we want to regularize equation (^) for the electronic transition dipole M. This equation 
is highly nonlinear: in addition to the dependence on M in the wave function, it also involves 
second and third powers of M. Once again, we may choose at which point to make the lineariza- 
tion. We prefer to linearize at an early stage, in order to obtain an equation of a simple form. 
Thus we define 



FM{x,t) = M\x)F{x,t) 

and write the functional 

ft 



GMix,t) = M^{x)G{x,t), 



(23) 



Jm 



dt' 







dxMix)GMix,t') + j dxFM{x,t') - QMit') 



+ aM / dxM'^ix). (24) 



Here, guit)) is the right-hand side of equation ^. The regularized solution of equation (g) is 
achieved by minimizing this functional. And we choose the linearization procedure by taking 
the variation of Jm in equation (|23) only with respect to the explicit dependence on M (that 
is, we treat Fm and Gm as independent of M). Then we obtain 



where 



dx'^M{x,x' ,t)M{x') + aMM{x) = 0M(x,t), 



^M{x,x',t)= / dt'GM{x,t')GA4{x',t'), 



QMix,t) = f dt'GMix,t') guit') - f 



gM{t')- / dx'FMix',t') 



(25) 

(26) 

(27) 



Finally, the integral equations (12 0| ) and (25) for the potentials and dipole moments and the 
Schrodinger equations (|^ and (||) for the components of the wave function form the full set of 
coupled equations for the unknown functions, with (/(t)), eo(*)) and ei(t) as input data. 

We conclude the presentation of the formalism with a schematic outline of the inversion 
algorithm which will be numerically implemented in a forthcoming work. First, the algorithm 
will start with trial functions for the potentials and dipoles to propagate the wave function from 
t = to t = At. Using trial functions at the first step is not an excessive demand for two 
reasons: (i) For a sufficiently small time increment At, the evolution of the wave function is 
mainly affected by the values of the potentials and dipoles in the region where uo{x) is localized, 
i.e., in the harmonic region of the ground potential surface; such information is usually known 
with reasonable accuracy, (ii) As more data becomes available, the initial trial functions will 
be replaced by those which match the measured fluorescence intensity. The second step will 
use the measured fluorescence intensity, the wave function components u{At) and v{At), and 
the initial trial functions for solving equations ( p5[ ) and (^) to obtain the next evaluation of 
the potentials and dipoles. These functions will be once again substituted into the Schrodinger 
equation to propagate the wave packet from t = At to t = 2At. The procedure will be repeated 
many times with new laboratory data incorporated at each time step. The recorded fluorescence 
intensity, {I{t)), contains information about the potentials and dipoles in the region where the 
wave packet is localized at moment t as well as where it was prior to that time. The sequential 
marching forward in time over the data track acts to refine the potentials and dipoles at each 
time step. 

4 Discussion 

This letter sets forth the formalism of a novel comprehensive approach to the inversion of 
molecular dynamics from time-dependent laboratory data. One of the main features of the 
proposed inversion method is initiation in the well-known ground state and use of external 
driving fields to excite the wave packet and guide its motion on the ground and excited potential 
surfaces. Different driving fields will induce different dynamics and may be more or less helpful 
for the inversion procedure. Consequently, we are left with the attractive prospect of choosing 
the driving fields to be optimally suited for assisting the extraction of unknown potentials and 
dipoles from laboratory data. This choice may be facilitated by a closed learning loop |[l^ in the 
laboratory, starting with a number of different trial fields. According to the inversion objectives, 
a learning algorithm will determine the best candidates and direct the fields to shapes which 
are best suited to produce these objectives. Natural objectives are to maximize the spatial 
region where the potentials and dipoles are reliably extracted. Physical intuition suggests that 
one may learn more about the Hamiltonian at a specific spatial point if the wave packet is not 
spread over the whole potential surface but is essentially localized in a narrow region around this 
point. Consequently, the driving fields best suited for the inversion will control the dispersion of 
the wave packet and guide its motion in a desired large spatial region. A numerical simulation 
of the algorithm, including its closed-loop learning features, will be the next step towards its 
laboratory implementation. 
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